Two-loop Critical Fluctuation-Dissipation Ratio 
for the Relaxational Dynamics of the 0{N) 
Landau- Ginzburg Hamiltonian 



Pasquale Calabrese and Andrea Gambassi 
Scuola Normale Superiore and INFN, Piazza dei Cavalieri 7, 1-56126 Pisa, Italy. 
e-mail: calabres@df.unipi.it, andrea.gambassi@sns.it 
(February 1, 2008) 



Abstract 



The off-equilibrium purely dissipative dynamics (Model A) of the 0{N) vec- 
tor model is considered at criticality inane = 4 — d>0 expansion up to 
O(e^). The scaling behavior of two-time response and correlation functions at 
zero momentum, the associated universal scaling functions and the nontrivial 
limit of the fluctuation-dissipation ratio are determined in the aging regime. 
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I. INTRODUCTION 



In recent years many efforts have been made in order to understand the off-equihbrium 
aspects of the dynamics of statistical systems. A variety of novel dynamical behaviors 
emerges when some kind of randomness is present in the system. Among them, one of the 
most striking is that of aging (see Ref. [1] and references therein). It has been pointed out [2] 
that they could also emerge in nondisordered systems if slow-relaxing modes are present. 
This naturally happens when the system undergoes a second-order phase transition at some 
critical temperature Tc- Indeed, consider a ferromagnetic model in a disordered state and 
quench it to a given temperature T > Tc [3] at time t = 0. During the relaxation a small 
external field h is applied at x = after a waiting time s. At time t, the order parameter 
response to h is given by the response function i?x(t,-s) = 6{(f):x,{t))/6h{s), where (p is the 
order parameter and (•) stands for the mean over the stochastic dynamics. Correlations of 
order parameter fiuctuations are interesting dynamical quantities as well. In the following 
we will focus on the two-time one, given by Cx{t,s) = (0x(^)0o(s))- The time evolution 
of the model we are considering is characterized by two different regimes: a transient one 
with off-equilibrium evolution, for t < tR, and a stationary equilibrium evolution for t > tn, 
where tR is the relaxation time. In the former a dependence of the system behavior on initial 
condition is expected, while in the latter time homogeneity and time reversal symmetry (at 
least in the absence of external fields) are recovered; as a consequence we expect that for 
tR < s,t, R^{t,s) = R^^{t - s), C^{t,s) = C^%t - s) where R"^ and are determined 
by the "equilibrium" dynamics of the system, with a characteristic time scale diverging at 
criticality (critical slowing down) . Moreover the fluctuation-dissipation theorem states that 

= -'J^. (1.1) 

If the system does not reach the equilibrium, all the previous functions will depend both 
on s (the "age" of the system) and t. This behavior is usually referred to as aging and was 
flrst predicted for spin glass systems [1,4]. The fluctuation-dissipation ratio (FDR) [2,4]: 

is usually introduced to measure the distance from the equilibrium of an aging system 
evolving at a fixed temperature T. A non trivial value for this ratio is also experimentally 
observed in some glassy systems [5]. 

In recent years several works [1,2,4,6-14] have been devoted to the study of the FDR 
for systems exhibiting domain growth [15], and for aging systems such as glasses and spin 
glasses. Xx(t, s) turns out to be a nontrivial function of t and s, in the low-temperature 
phase of all these systems. In particular, analytical and numerical studies indicate that the 
limit 

X~ = lim lim Xx=o(t, s), (1.3) 

s— >oot— >oo 

vanishes throughout the low-temperature phase both for spin glasses and simple ferromag- 
netic systems [8-10,12,13]. 
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Only recently [2,16-21] attention has been paid to the FDR, for nonequilibrium, nondis- 
ordered, and unfrustrated systems at criticality. From general scaling arguments one would 
expect that the critical response function scales as [19-23] : 

i?x=o(i, s) = Anit - sr-^/%t/sfj^R{s/t) , (1.4) 

where a — {2 — r] — z)/z and 9 is the initial-slip exponent of the response function, related to 
the initial-slip exponent of the magnetization 9' and to the autocorrelation exponent Ac [24] 
by the relation [22] 

0' ^0 + z-\2-z-r))^z-\d-Xc). (1.5) 

In recent works the notion of local scale invariance has been introduced as an extension 
of anisotropic or dynamical scaling (see [25] and references therein). Assuming that the 
response function transforms covariantly under the constructed group of local transforma- 
tions, it has been argued [26] that J^R{s/t) — 1. Under the same assumption, the full spatial 
dependence has been also predicted [25] 

R^{t, s) = R^=o{t, s)<f (|x|/(i - sy/'), (1.6) 

where $(-u) is a function whose convergent series expansion is explicitly known [25]. For the 
correlation function and its derivative no analogous prediction exists. One can only expects 
from general Renormalization Group (RG) arguments that [19-23] 

Cx=o(i, s) = Ac{t - sY+'-''^{t/sf-'Tc{s/t) , (1.7) 
dsC^=o{t, s) = Aacit - sY-^l^t/sYTacisIt) , (1.8) 

with the same 9 and a as in Eq. (1.4). The functions Tc{v), ^dc{v) and J-r{v) are universal 
and defined in such a way ^c(O) = ^9c(0) = •^ii(O) = 1, while the constants Ac, Adc and 
Ar arc nonuniversal amplitudes. 

From the scaling laws of above, it has been argued that X°° is a new universal quantity 
associated with the given nonequilibrium dynamics [19,17,20], and, as such, it should attract 
the same interest as critical exponents. Given this universality, it is worthwhile to compute 
X°° for those mesoscopic models of dynamics which have the same critical behavior of some 
lattice models considered so far in the literature. 

Correlation and response functions were exactly computed in the simple cases of a Ran- 
dom Walk, a free Gaussian field, and a two-dimensional XY model at zero temperature and 
the value X°° = 1/2 was found [2]. The analysis of the d-dimensional spherical model gave 
X°° ^ 1-2/d [17], while = 1/2 for the one-dimensional Ising-Glauber chain [16,19]. 
Monte Carlo simulations have been done for the two- and three-dimensional Ising model [17], 
finding = 0.26(1) and X°° ~ 0.40 respectively. The effect of long-range correlations in 
the initial configuration has been also analyzed for the rf-dimensional spherical model [27]. 

Only in a recent work [21] field-theoretical methods have been applied to determine the 
FDR and the scaling forms of the response and correlation functions up to the first order in an 
e expansion, for the purely relaxational dynamics of the 0{N) Ginzburg-Landau model. This 
field-theoretical model has the same symmetries (and thus the same universal properties) as 
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a wide class of spin systems on the lattice with short range interactions (see [30], or [31] for 
a recent review). 

In [21] the following quantity, related to the FDR, was introduced in momentum space 

where i?q(t, s) and Cq(t, s) are the Fourier transforms (with respect to x) of Rx{t,s) and 
Cx(t, s) respectively. It was argued that the zero-momentum limit 

X^o = lim lim Xa=o{t,s) (1.10) 

is equal to the same hmit of the FDR (1.3) for x = 0, i.e. = to all orders [21]. 

This fact allows an easier perturbative computation in momentum space of the new universal 
quantity X°°. 

The extension of these results to two-loop order is very important not only from a 
quantitative point of view. In fact, in the past, when new scaling relations have been 
proposed, several times they resisted to the test of the first order in the e expansion, but 
not to higher order calculations. Classical examples may be found in the context of surface 
criticality (see e.g. the review [28], p. 116, and references therein) and in the case of 
anisotropic scaling at Lifshitz points [29]. This is a further reason to present here the second 
order computation of the scaling form for the zero-momentum response function and the 
FDR for the purely dissipative relaxation of the 0{N) model. 

The paper is organized as follows. In Sec. II we briefly introduce the model. In Sec. Ill 
we evaluate the zero- momentum response function Rq=o{t, s) and in particular we derive its 
scaling form. In Sec. IV we compute the FDR up to the second order in e and we derive a 
scaling form for dsCci=o{t, s). In Sec. V we summarize and comment our results and discuss 
some points that need further investigation. In the Appendices A and B we give all the 
details to compute the zero-momentum Feynman integrals. 



II. THE MODEL 



The time evolution of an iV-component field (^(x, t) under a purely dissipative dynam- 
ics (Model A of Ref. [32]) is described by the stochastic Langevin equation 



a,(^(x,t) = -n-^^ + e(x,t). 



(2.1) 



Mx, t) 

where Q is the kinetic coefficient, 'C(X)^) ^ zero-mean stochastic Gaussian noise with 

(e,(x,i)e,(x',0) = 20(5(x-x')5(i-0<^i,-, (2.2) 

and Tl[(fi\ is the static Hamiltonian. It may be assumed, near the critical point, of the 
Landau-Ginzburg form 



(2.3) 
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Instead of solving the Langevin equation for ^p{^) and then averaging over the noise 
distribution, the equihbrium correlation and response functions can be directly obtained by 
means of the field-theoretical action [30,33] 



(2.4) 



Here (^(x, t) is an auxiliary field, conjugated to the external field h in such a way that 
7i[(/9, h] = Tifv?] — / d'^xhip. As a consequence, the linear response to the field /i of a generic 
observable O is given by 

^^^^ -n{^(^,s)0), (2.5) 



5/i(x, s) 



for this reason (^(x, i) is termed response field. 

The effect of a macroscopic initial condition v'o(x) = V'(x, t = 0) may be taken into 
account by averaging over the initial configuration with a weight e~^°^'^°^ where, for example, 

i/oM = /d'^x^((^o(x)-a(x))^ (2.6) 

that specifies an initial state a(x) with gaussian short-range correlations proportional to Tq"^. 
Any addition of anharmonic terms in Ho[ipo\ is not expected to be relevant as long as the 
harmonic term is there (as in the case when the initial state is in the high temperature phase). 
Instead, an initial condition with long-range correlations may lead to a different universality 
class, as e.g. shown for the d-dimensional spherical model with non-conservative dynamics 
[27]. 

Following standard methods [30,33], the response and correlation functions may be ob- 
tained by a perturbative expansion of the functional weight e~^'^['^''^l+^°['^°l^ in terms of the 
coupling constant go (appearing in the vertex gQ(fi^(f/3\). The propagators (Gaussian two- 
point functions of the fields (p and (p) of the resulting theory are [22] 

((^,(q, s)ip,{-q, t))o = 5,,Rl{t, s) = 5,, e{t - s)G'(t - s), (2.7) 



{ipi{q, s)^j{-q, t))o = SijC^{t, s) 
where 



+ ro 



^-.|)+(^^-l) G{t + s) 



(2.8) 



G{t) = e-"(«'+''°)*. (2.9) 



The response function Eq. (2.7) is the same as in equihbrium. Eq. (2.8), instead, reduces to 
the equilibrium form when both times t and s go to infinity and t — t — s is kept fixed. In 
the following we will assume the Ito prescription (see [23,34], [30] and references therein) to 
deal with the ambiguities that arise in formal manipulations of stochastic equations. Con- 
sequently, all the diagrams with self-loops of response propagator has to be omitted in the 
computation. This ensures that causality holds in the perturbative expansion [22,23,33]. 
Prom the technical point of view, the breaking of time homogeneity makes the renormaliza- 
tion procedure in terms of one-particle irreducible correlation functions less straightforward 
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FIG. 1. Two- loop Fcynman diagrams contributing to the response function. Response propaga- 
tors are drawn as wavy-normal lines, whereas correlators are normal lines. A wavy line is attached 
to the response field and a normal one to the order parameter. 

than in standard cases (see Ref. [28,22,23]). Thus the computations will be done in terms 
of connected functions. 

From the expressions above, we can compute the FDR for the Gaussian model [2,21]: 

X^it, s) = (^1^ j ' = (l + e-^^(^^+'-")^ + nq'T^'e-'''^'^'^''-^>y' . (2.10) 

When the model is not at its critical point, i.e. tq (x T — Tc ^ 0, the limit of this ratio 
for s — > oo is 1 for all values of q, according to the idea that in the high-temperature 
phase all modes have a finite equilibration time. In this case equilibrium is approached 
exponentially fast in time and as a consequence the fluctuation-dissipation theorem applies. 
For the critical model, i.e. Tq = 0, if g 7^ then the limit ratio is again equal to one, whereas 
for g = we have A'°^Q(t,s) = 1/2. We can argue that, in the Gaussian model, the only 
mode characterized by aging, i.e. that "does not relax" to the equilibrium, is the zero mode 
in the critical limit. 



III. TWO-LOOP RESPONSE FUNCTION 

In this Section we compute, up to the second order in a loop expansion, the critical 
nonequilibrium response function at zero external momentum for the model described in 
the previous Section. We use here the method of renormahzed field theory in dimensional 
regularization with minimal subtraction of dimensional poles. Up to the second order in 
perturbation theory there are four connected Feynman diagrams (without self-loops of re- 
sponse propagator) that contribute to the response function. They are depicted in Figure 1. 
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In terms of these diagrams and as a function of the bare coupUngs and fields (denoted in 
the following with 9?^, (ps)-, the zero-momentum bare response function s) is given by 
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+5o 



iV + 2 
6 



(iV + 2)- 
18 
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2,3 



+ 0{gl). 



(3.1) 



In the following we assume t > s for simplicity. We also fix Tq ^ 



iQ — 0, since Tq"^ is an 
irrelevant variable (in the RG sense) and thus it affects only the corrections to the leading 
scaling behavior [22,23]. Using the results reported in the Appendix A, we get 



RB{t,s) = 1+ ^0 



N- 



24 



^°^7 + 2 



(7e + log 2 + log t) log - - i log^ - 

s 2 s. 



+ 



144 



-9o 24 



(iV + 2)2 t 



log' - + {N + 2y 



Q + log2 + 7b + logi) log ^ + ^ log' ^ I 



- (^log - + log - j + log - + logt + 7eJ - - log - 

+ 0{go,go6,goe^) , 



+(log(t-.) + 7i.)log^-^ 



(3.2) 



where ^0 = ^dgo, — 2/((47r)'^/'r(d/2)) and f{v) is a regular function defined in Eq. 
(A28). To lighten the notations we set Q = 1 in the previous equations. The dependence on 
Q of final formulas may be simply obtained hj t Qt, where t is the generic time variable. 

In order to cancel out the dimensional poles appearing in this function, we have to 
renormalize the coupling constant according to [30] 



^0 



1 + 



N + 8g 
6 7 



~9 + 0{g% 



(3.3) 



and the fields (p and ip via the relations [33] ipB — Zl/'^ip, (pB — Z^'^ip, so that 



R{t,s) = {Z^Z^)-^'^RB{t,s) 



RB{t,s). 



After this renormalization, R{t, s) is a regular function of dimensionality also for e 
The critical response function is now obtained by fixing g at its fixed point value 



6e 



N + 8 



3(3Ar + 14) 

1+ ^TT TTTT^e 



leading to 

R{t,s) = 



N + 2 ^ t e' 
' + ^4(]vT^^°S;+4 



(AT + 8)2 
6(A^ + 2) fN + 3 



(3.4) 
^ 0. 

(3.5) 



(A^ + 8)2 \N + S 
3{N + 2 



+ log 2 ) log - + , j,, log' - 

s 



6(A^ + 2), 4, , 

■(TVW '3^°'^'-'^ + 2(]vT8p 



/(s/t)-47Elog 



s 8(Ar + 8)2 

4^ 



(3.6) 
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FIG. 2. Plot of the two-loop contribution to the universal functions Fji{v) (see Eq. (3.11)) and 
Fqc{v) (see Eq. (4.9)). 



Note that the nonscaling terms, like logtlogt/s (appearing, for example, in i? ' , see 
Eq. (A27)), cancel each other out when the coupling constant is set equal to its fixed point 
value. Eq. (3.6) agrees with the expected scaling form in momentum space (analogous to 
that in real space, Eq. (1.4)) 



R{t,s)^AR{t-sY{t/sfFR{s/t), 
with the well-known exponents [22,23,30] 



a 



N + 2 e 
iV + 84 
2-77- 
z 



1 + 



6e 



N + 8 \N + 



N + 3 



+ log2 



3(A^ + 2) , 4 2 3, 



and the nonuniversal amplitude 

2 3(iV + 2) 



1 + 



8{N + 8)^ 

For the new universal function Fji{v) we find 

2 3(iV + 2) 



/(0)-47^1og^)+O(6=^) 



Fr{v) = 1 + 



8{N 



{f(v)-m) + 0{e'), 



(3.7) 

(3.8) 
(3.9) 



(3.10) 



(3.11) 



A plot of the quantity f{v) — /(O) (defined in the Appendix A Eq. (A28)), that completely 
characterizes the out-of-equilibrium corrections to the mean-field behavior up to the second 
order in the e-expansion, is reported in Fig. 2. Due to the small prefactor (e^/72 for the 
Ising model, N — 1), it might be very hard to detect these corrections in numerical and 
experimental works, as it happens for the corrections to the mean-field behavior of the 
static [31] and equilibrium dynamics [35] two-point functions. 
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FIG. 3. Diagrammatic trick. 



IV. TWO-LOOP FLUCTUATION-DISSIPATION RATIO 

In this Section we evaluate the FDR up to the order e^. We do not compute the full 
two-point correlation function C{t, s), since only dsC{t, s) is required to determine the FDR. 
This derivative may be computed by using the following diagrammatic identity. 

Each amputated diagram Di{t, s) (with label i) contributing to the response function, 
also contributes to the correlation one in two diagrams, as graphically illustrated in Fig. 3. 
Taking into account the explicit form of the propagators (sec Eqs. (2.7) and (2.8)) for = 
and causality (which also implies that Di{t,s) oc 9{t — s) apart from contact terms) it is 
easy to find that 

roo 

dsCi{t,s)^2Ri{t,s) + 2 dt't'Di{t',s), (4.1) 

Jo 

where Cj(t, s) is the contribution of this diagram to the correlation function, Ri{t,s) the 
contribution to the response function, and Di{t', s) the common amputated part. 




s t' t" t 

^2,5 

FIG. 4. Diagrams contributing only to the correlation function. 
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Relation (4.1) is nothing but a particular case of a relation following an algebraic identity 
for the functional integral, i. e. 

= 1 [d(/^d(^d(^o] , J_^ |<^(x', i)e-^[^'^]-^o[^ol I , (4.2) 



5ip{x,s) I- J 

with t > s > 0. At criticality (i.e. ro = 0, using dimensional regularization) we get in 
momentum space 



9o 



{ds - q )(<^(-q, t)<^(q, s)) = 2(<^(-q, t)(^(q, s)) - ■^('^(-Q, ^)<^ (Q, «)) , 



(4.3) 



which, in the limit q^ 0, is diagrammatically expressed by Eq. (4.1) as far as common 
amputated contributions to response and correlation functions are concerned. 

Diagrams contributing to the correlation function, but not to the response one do exist. 
They have to be computed without taking advantage of this identity. At two-loop order 
there are two of them, as in Fig. 4. 

Summing the six contributions to the correlation function we finally arrive to the expres- 
sion: 



dsCB{t,s) 



= R{t, s) - go 



6 



V 6 



{dC) 



6 



1 

+ 2 



N + 2 
18 



N ■ 



6 



{dC) 



2,5 



+ 0{gl) . (4.4) 



Considering the explicit expression for the diagrams given in the Appendix B one obtains 
the derivative of the bare correlation function. This bare quantity is renormalized using 
equations (3.3), (3.4) and 

n = ZaftB with Zn = ( ^ 1 , (4.5) 

so that, taking into account the ^2 wc set equal to 1 in the previous relations, 

dsC{t, s) = ZnZ-' dsCsit, s) = (Z^Z^)-'/^ dsCB{t, s) . (4.6) 

The expression of dsC{t, s) in terms of the renormalized coupling has a multiplicative redef- 
inition of its amplitude at the first order in g. Considering the fixed point value for g (cf. 
Eq. (3.5)) one finally obtains 



1 + e 



X |l + e 

{N + 2f 
8{N + 8y 

X 1 1 + 



N + 2 
4{N + 8) 
N + 2 



+ e- 



3(iV + 2)(3iV+ 14) 
4(iV + 8)3 



- - 
4(A^ + 8) s ^ 4 



6{N + 2) fN + 3 



{N + sy \N + S 



+ log 2) log ^+ 



s 

N + 2 



t 6{N + 2) 4 



3, 4, t-s 3, 4 , ^ 3 , , A^ + 2 

,4 -3''^ ITS- '^2 3 + ^^^/^^ + 8^^^/^) + 



" t) 2 - ^ (1 + log + l V I + 1^,(1/ A)] I + 0{e') , (4.7) 
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where the function f{v) and $(f) are defined in Eqs. (A28) and (B12) respectively, and 
Li2 is the dilogarithmic function whose standard definition is recalled in Eq. (A15). Note 
that also for dsC{t, s) all the nonscahng terms cancel out when the couphng constant is set 
equal to its fixed point value. This result agrees with the scaling form in momentum space 
(analogous to Eq. (1.8)) 



dsC{t, s) = Aocit - sr{t/sYFac{s/t), 



(4. 



with the same a and 9 as given in the previous section and a new universal scaling function 
Fqc{v) given by 



Fac{v) 



, 3(iV + 2) 
8(A^ + 8)2 



2 log 1 log i-p^ + ^mv) - m) + fiv) - /(O) 



O(e^) 



(4.9) 



A plot of the loop corrections in the above expression (apart from the factor gl^l^^^^2 appearing 
also in -Fr(v)) is shown in Fig. 2. As already noticed for -Fr(v), effective corrections to mean- 
field behavior are quantitatively very small for Fqc{v)- 

Taking the long time limit (according to Eq. (1.10)) of both the correlation and response 
functions one obtains the limit of the critical fiuctuation-dissipation ratio we are interested 



m: 



oo \ — 1 
q=0J 



= 1 + 



-e + e 



N 



4(A^ + 8) (AT + 8)2 



A^ + 2 3(3Ar+14) 



8 



4(A^ + 8) 



+ 0(e3), (4.10) 



with 



3 3 23 
c = - ^ + ^ log 2 (2 + 1 1 log 2 - 3 log 3) - - log^ 3 

+ ^Li2(l/4) - |Li2(l/3) + yLi2(3/4) - ^Li2(8/9) 



-0.0415... . 



(4.11) 



We note that the contribution of c to the FDR is quite small. For example, with A^ = 1 the 
sum of the first two terms in brackets is ~ 1.8, which is about 45 times larger than c. 



V. CONCLUSIONS AND DISCUSSIONS 

In this work we studied the off-equilibrium properties of the purely dissipative relax- 
ational dynamics of an A^- vector model in the framework of field theoretical e-expansion. The 
results presented here extend those of Ref. [21]. The scaling forms for the zero-momentum 
response function and for the derivative with respect to the waiting time of the two-time 
correlation function read 

R{t, s) = Anit - snt/sfFnis/t) , (5.1) 
d,C{t, s) = Agcit - sy{t/sfF9c{s/t). (5.2) 
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0.5 1 1.5 2 2.5 e 



d 4 3 2 

FIG. 5. X°° as a function of the dimensionality d = 4 — e for several N. For each N the upper 
curve is the [2,0] Padc approximant and the lower one the [0,2]. The exact result for A'^ = oo is 
reported as a solid line. The numerical Monte Carlo values for the Ising Model in two and three 
dimensions are also indicated (for the latter, there is no indication about the error). 

The universal functions Fji{s/t) and FQc{s/t) are given in Eqs. (3.11) and (4.9) respectively. 
In both cases the corrections to the Gaussian value 1 is of order e^. In principle these cor- 
rections should be detectable in computer and experimental works, but being quantitatively 
very small, they are hard to observe. We would remark that this fact does not mean that 
aging effects in these models are weak compared with the analogous phenomena in glassy 
systems. In fact aging manifests itself in the full scaling forms (e.g. ^ 7^ 0) and in the 
violation of fluctuation-dissipation theorem, i.e. in 7^ 1 in a quantitative way. 

We note that the i?(t, s) we found agrees with the general RG form, but at first sight 
it is not compatible with the Fourier transform of Eq. (1.6). This naive comparison should 
be done very carefully because it involves a Fourier integral which could be divergent. The 
analysis of the full g-dependence of -Rq(t, s) may give some insight into this problem. This 
dependence has been already carried out up to 0{e) [21], but it is very hard to determine 
it up to two loops. In other dynamical universality classes this discrepancy already arises 
at 0(e). The computation of the full g-dependence in these cases seems to be simpler and 
may provide some useful hints [36]. 

We computed the FDR A'q=o for general iV, cf. Eq. (4.10). As shown in Ref. [21] this 
quantity for zero momentum has the same long-time limit as the standard FDR X^. Using 
this equality we may compare our result with those presented in the literature. 

In the limit 00, Eq. (4.10) reduces to = l/2-e/8-eV32-FO(e^), in agreement 
with the exact result for the spherical model = 1 — 2/d [17]. 

The formula for general (cf. Eq. (4.10)) allows us to make quantitative prediction for 
a large class of systems. In Fig. 5 we report the dependence of on the dimensionality 
at fixed N, while in Fig. 6 we show the dependence on N at fixed d ~ A — e. For each 
model we report two values: one is obtained by direct summation (Fade approximant [2, 0]) 
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0.5 




d = 2 



2 4 6 8 10 N 

FIG. 6. N dependence of for d = 2, 3. The upper curve is the [2, 0] Padc approximant and 
the lower one the [0, 2]. The dotted Une is the exact result for = oo in o? = 3 {X°° = 1/3) 

and the other by "inverse" summation (Fade approximant [0, 2]). We do not show the [1, 1] 
approximant, since it has a pole in the range of e we are interested in. Prom these figures 
some general trends may be understood: 

• decreasing the dimensionality, X°° always decreases, at least up to e = 2 (for the 
one-dimensional Ising model the value X°° = 1/2 is expected [17]); 

• increasing N, X°° decreases, approaching in a quite fast way the exact result for the 
spherical model; 

• for = oo the curve of the [0, 2] approximant reproduce better the exact result in 
any dimension with respect to the [2, 0] approximant. 

The last point suggests the use of the [0,2] value as estimate of X°°, also for physical 
A^. We quote as indicative error the difference between the two approximants. Using this 
procedure, we obtain X°° — 0.429(6) for the three-dimensional = 1 model, compared 
to ~ 0.46 found at one- loop [21], in very good agreement with the Monte Carlo simulation 
value ~ 0.40 for the three-dimensional Ising Model [17] with nonconservative (heat-bath 
Glauber) dynamics. Considering e = 2 one obtains X°° = 0.30(5) for N = 1, improving the 
one loop estimate ~ 0.42 in the right direction towards the Monte Carlo result = 0.26(1) 
for the two-dimensional Ising Model with Glauber dynamics [17]. 

Using our results we can give predictions of X°° for systems that have not yet been ana- 
lyzed by numerical simulations. We estimate X°° = 0.416(8) for the three-dimensional XY 
model and X°° = 0.405(10) for the three-dimensional Heisenberg model. These predictions 
may be tested by numerical simulations extending the results quoted in [17]. 

There are also several open questions that need further investigations. For example a 
"rigorous" proof of the fact that the FDR is exactly 1 for all modes with (somehow 
related to the presence of a mass gap) has not yet been given. Then one might ask how these 
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theoretical results (scaling forms, relaxing modes etc.) change if one changes mesoscopic 
dynamics (e.g. with conserved quantities), or when more complex static Hamiltonians are 
considered, e.g. those with disorder or frustration, or when different initial conditions (e.g. 
with long-range correlations) are considered. We will consider in forthcoming works the 
0{N) model with Model B and C dynamics [36] and the purely dissipative relaxation of the 
Ising model with quenched random impurities [37]. 
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APPENDIX A: CONNECTED DIAGRAMS FOR THE RESPONSE FUNCTION 



The four diagrams contributing to the response function up to the two-loop order are 
reported in Fig. 1. The one- loop diagram was already discussed in Ref. [21]. The expression 
obtained there for the critical bubble (i.e. for the IPI part of the diagram) is 



Beit) 



C%t) 



(2t) 1-^/2 V{d/2-l) 



{27rY ^'''"'"^ d/2-1 (47r)'^/2 
Thus the full connected 1-loop diagram for the response function is given by 



(Al) 



1 



^'2^^/2(2-^/2)^^' "' ^ 
t e . _ . . . t 1 . o t 



log - + - (7e -f- log 2 log t) log - - - log - 
s 2 V s 2 s 



+ 0(e' 



(A2) 



From these one-loop expressions, it is quite simple to compute the two-loop integrals 
R"^'^ and i?^'^ of Fig. 1. Indeed the two-loop critical bubble (the IPI part of R"^'^) can be 
computed in terms of Bc{t) as 



b{d) 



(A3) 



where 



b{d) 



r^{d/2-i) 

2d-i 



{i-d)T\2-d/2) 
2r(4 - d) 



-- [l + e(log2 + 7s) + 0(e2)]. (A4) 



By means of this expression, we compute R^'^ taking into account the external legs with 
q = 



(4-d)2 ' 

that near four dimensions has the following series expansion 



(A5) 
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t (I 



s \e 



log - - + logi + log2 + 7e + - log^ - 



t 



0{e). 



(A6) 



The computation of R ' is simple once the expressions for R ' and Bc{t) are known. 



Indeed, from Eq. (A2), it is obtained: 

R^'\s,t) = J^dt'R^'\t,t')B,{t') = TV, 



^2-d/2 _ g2-d/2 
4-d 



(AT) 



that is, expanding in e. 



t2 1 i„„2 t 



R'''is,t)^N!-\og'- + 0(e). 
32 s 



(A8) 



The last diagram R^'^ is more difficult to be worked out and it requires a long calculation 
which main steps are described in what follows. First of all we evaluate its IPX contribution 
called Oi{t, s) 



Oi{t,s) 



d% f d'^Ct ^0 x^O u .\ DO 



C'Jt,s)Cl{t,s)Rl^Jt,s) 



^ {27iY {27iY qjqr ' 



e{t - s) 



{t-sf-''Jd{lA) + {t + sf-''Ja{l, 



2-d 



t - S 
t + S 



-2{t-sY-''Jd 



t + s 



, (A9) 



with 

Jdia,b) 

and 

Fdix) -- 
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(27r)<^i {27rY 



{b+l){a + b)J ' 
(AlO) 



Tid/2) 



dss'^/2-2gW4r(0, s) = 



T{d/2)T{d/2 - 1) 
2{d - 2) 



d d 



d X 



2' 4 



In particular, for our calculations, we are interested in the limits 
log(l - x/4) 

-1 



X 



F4_,(l)=log- + e 



1 



7e - - log - - - lo; 



3 4^32 



Fd{0) = 



T\d/2 - 1) 



(All) 

(A12) 
(A13) 
(A14) 



Here Li2{z) is the standard dilogarithm, defined as 
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Li2(^) = E 



k=l 



A;2 ■ 



(A15) 



The final expression for Oi{t, s), in generic dimension, is 



Oi{t,s) 



2d-2 



t - s 



-2{t{t-s)Y-'"'Fa 



t - s 



Tlie full connected diagram is thus given by the following expression 

R^'^s, t) = f dt" /* dt'0,{t', t") = Nl{A,{s, t) + A^is, t) - 2A^{s, t)), 

Js Jt" 



(A16) 



(A17) 



where 



A,{s,t)^2'-''Fail) 



(t-s) 



A-d 



(A18) 
(A19) 



(4-rf)(3-rf) ' 

Js/t Jy 

As{s, t) = 2''-H^-'' r dy y''-'' C dz /"^(l - zf-'^'^Fdil - z) = 2'-H^-''h{s/t) . (A20) 

Js/t Jy 



The evaluation of the two functions Ii{v) and l2{v) is rather cumbersome but algebraically 
trivial. After some calculations one gets 



I^^y) = £:(!!Z^J)[iog^(81og2 - 61og3) + Mv) + 0(e)], 



4 



r\d/2-i) 



--logv- log^ t; - (1 - 6 log 2 + 3 log 3) log^; + f2{v) + 0(e) 
where fi{v) are given by 

/i(^;)/4 = logt; r dzF^{{l-z)^)+ [' dz log z F4{1 - zf), 

Jo Jv 



h{v)/A = logv / dz 



F^{l-z)-l , /•! logz 
' ' dz 



l-z 



l-z 



[F4(l - ^) - 1] + / dz 



1 log(l-^) 



and in particular these are regular functions in the limit v — > 

/i(0) = V 2 + V \ + 3Li2(l/4) - 4Li2(2/3) , 
/2(0) = -^ + ^V^-Li,(l/4). 

Inserting all these contributions in Eq. (A17), we get 

-- (^log - + log ^ ) - log -(log(i - s) + 7fi) 



(A21) 
(A22) 

(A23) 
(A24) 

(A25) 
(A26) 



3' 
t 1 
~s ■ 2 



-(2 + 7^; + logi) log - + - log2 - + 



+ 0(e). 



(A27) 



16 



with 



fiv) = - 2f2{v) - log ^(2 + log 12) - 2Li2(l/4), (A28) 

/(O) = y - 2 log ^ + 3 log2 2 - V I + 3Li2(l/4) - 4Li2(2/3) 

= 0.663707.... (A29) 



APPENDIX B: CONNECTED DIAGRAMS FOR THE FDR 



In this Appendix we evaluate the rest of the diagrams required for the computation of the 
FDR. Wc do not evaluate the full integral for the correlation function, since we make use of 
the trick explained in details in Section IV. For this reason we consider first those diagrams 
contributing also to the response function and we evaluate only their extra-contributions 
(given by dt"t"Di{t", s) in Eq. (4.1) and denoted with the subscript "e" in what follows) 
to the derivative of the correlation function. For the first three diagrams these contributions 
are very simple: 



-\ - ^(logs + 7B + log2) 



(dcr/ = i; mnB^is) = = ^ + o(e) 

{dCff = sB,^{s) = -Nl\{- + + log2 + logs) + 0(e) . 

o e 



(Bl) 
(B2) 
(B3) 



The fourth contribution is less simple 
Jo 



4-d 



FJl] 



(4-ci)(3-ci) 

+ dz zFd{{l - zf) -2 f^dz z{l - zy-'^/^Fdil - z) 
Jo Jo 



(B4) 



Using now the explicit form for ^4(0;) given in Eq. (A12), one obtains 



4{dC) 



2,3 

e 



-(^ + logs + 7B + ^) (log^ + 1) 



Li2(l/4) + log ^ (J log ^ - log 2 ) + 0(6) . (B5) 



The diagrams whose amputated part does not contribute also to the response function 
are shown in Fig. 4. The sunset-type diagram (dC)'^''^ is quite difficult, thus wc first compute 
its IPI part 02{t, s). Introducing Qs = qi + q2, this contribution is given by (for t > s) 



02{t,s) 



d'^Qi f d'^q2 
{2TrY 
d'^qi 



{2tt 
d'^q2 



A{t,s)Cl{t,s)Cl{t,s) 



{2nYJ (27r) 



n 



qf 

,3-d 



3-d 1 



,3-d i 



(B6) 
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with A = t — s, a = t + s, and 

1 f d% r d% 1 1 1 



J (2nYJ (2n) 



g-9?-9i-a:(9l+92)^ 



NiJ {27:YJ {27:Yqlql{qi + q2Y 
r(d/2 - l)r(d/2) ^oo f\^^d/2^2 



4 ^ (1 + u) 

In the following we are interested in the limits 

1 2{l+x) 1 1 + 2.T 



vu 



{l + uf 



l-d/2 



31og^Q + 7i?) + ^log'3 + 



Li2(8/9) 



Introducing these results in the expression for the connected diagram 

{dC)''' = dt'O^it', s) = O^is, t') + 02{t', s), 



one finds 



4-(i 
+3 Tdi' 

^0 



+3 



{s + t'f-''K 



't'-s' 
t' + s, 



it' 



\3-d 



's + t" 
t'-s, 



(25^^'^ +{t- s)*-" -{t + s^-^) 

(Jj 

+3{2sr-' \ f dy{l + vY-'iKaiy) - y'-^Kail/y)] 
Jo 

+ dy{l - yr-'[Kd{y) - y'-'K,{l/y)] 
Jo 

^ log ^ (; + logs + ^ log ^ + 7e) + + 0(e) , 



where 



$(^) = 2i^^(l)-^7^;log^ + 



+3 



lo 1 



^ (i^4(|/) - iK.il/y)) + 



1 -y V 



In particular we are interested in the limit v ^ 0, given by 
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$(0) = 2X1(1) + 6 - 
Jo 1 



dy 



KM - -K, ( - 

y \y. 



3 4 
27.log- 



3, 4 39, 2^ 



^ log 2 + log 3 - ^ log2 3 - §Li2(l/3) 



+ ^Li2(3/4) - ^Li2(8/9) = -0.24889 . . . 
8 8 

Now the only diagram left is C^'^ of Fig. 4. It is given by 

C^'\t,s) = j dt'dt''R%{tX)Bc{nC'^=oit'',tW^ 



r^{d/2 - 1) 



■2d-^{3-d/2){2-d/2) 
Its derivative with respect to s, near four dimension, is 



g2-d/2 

5-d 



J,-d/2 



1 



log - + 1 

s 



+ 0{e) 



(B13) 



(B14) 



(B15) 
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